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Abstract: We investigate the performance of bends and splitters in 
graphene nanoribbon waveguides. Although the graphene waveguides are 
lossy themselves, we show that bends and splitters do not induce any 
additional loss provided that the nanoribbon width is sub- wavelength. We 
use transmission line theory to qualitatively interpret the behavior observed 
in our simulation. Our results pave a promising way to realize ultra-compact 
devices operating in the terahertz region. 
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1. Introduction 

Graphene has attracted considerable attention due to its unique electronic and optical proper- 
ties (U m m m O . The high carrier mobility approaching 200,000 cm^V~^s~^ offers the pos- 
sibility of building high-speed optoelectronic devices, such as photodetectors ||6l and optical 
modulators |7|. Despite being a single atomic layer, graphene can couple strongly to light, en- 
abling a constant absorbance of tt^z = 2.3% in the infrared-to-visible spectral range JS], where 
a = jlic is the fine-structure constant. More importantly, the optical properties of graphene 
can be substantially modified by changing the Fermi energy Ep (91 [T0| . For a slightly doped 
graphene {Ep < ;?z()l)/2), the conductivity of graphene is dominated by the interband contribu- 
tions. In the terahertz range when > fico/l, the intraband contribution becomes relevant, 
thus allowing graphene to support graphene-plasmon polaritons (GPPs), i.e., the collective os- 
cillation of carriers. Compared to conventional plasmonic materials, e.g. gold, doped-graphene 
plasmon has promising properties including tunability, extreme field confinement and low prop- 
agation loss. 

Graphene-plasmon polaritons supported by doped graphene have been widely investi- 
gated (HI [121 [IS [Ml [151 • It has been shown that remarkable optical absorption enhance- 
ment (T2I [131 [m in the graphene can be achieved by patterning the single sheet of doped 
graphene. Moreover, the GPP can also be confined laterally, e.g., in graphene ribbons. Re- 
cently, much attention has been focused on investigating propagation properties of GPP in 
graphene nanoribbons |[l71[l8l, which can be used as a promising candidate to construct future 



high-speed compact-devices working in the terahertz region. Obviously it is indispensable to 
have graphene nanoribbon bends and splitters as the basic structures in these compact devices. 
However, there is a lack of characterization of the GPP in these geometries. In this paper, we 
investigate the performance of bends and power splitters in graphene nanoribbon waveguides, 
particularly focusing on the question whether these bends and splitters will induce reflection or 
excess scattering loss on top of the propagation loss in the nanoribbon waveguides themselves. 
Our analysis in this paper relies on a classical electromagnetic description where graphene is 
described by a local conductivity, and the results presented here do not include the microscopic 
details of the edges of graphene nanoribbons [[TtI . 
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Fig. 1. (a) Effective refractive indices of GPP modes as a function of wavelength for free- 
standing graphene ribbons of 20 nm (the solid line) and 30 nm (the dashed line) widths. 
Inset illustrates a cross-section graphene nanoribbon waveguide with the width d, where the 
axis indicates light propagation along the z direction, (b) Propagation lengths as a function 
of wavelength for the freestanding graphene ribbons. Inset shows field distributions of the 
GPP mode at A = 15 jim for the case of d = 30 nm. 



2. Mode properties of GPP in a graphene nanoribbon 

Let us first investigate mode properties of GPP in a freestanding graphene nanoribbon, as 
shown in the inset of Fig[T] (a), with an ultra-thin width d<50 nm. The inset in Fig. (Ha) il- 
lustrates a cross-section of the nanoribbon and its corresponding axis. Due to the translational 
symmetry in the z direction, the electric solution of the nanoribbon waveguide has the form 
E(r, (O^t) = 'Ej{x^y)exp{ikz{(i))z — icot), where kz{co) is the wave vector in the propagation direc- 
tion z. For a waveguide in which the material is lossy, kz{co) is complex. An effective refractive 
index ^eff of the waveguide mode is defined by n^ff = Re{kz{co)) /ko, and the propagation length 
L of the mode is given by L = l/[2Im{kz{co))], where ko = 27l/^ (wavelength in vacuum Aq). 
In the simulations, the graphene ribbon is modeled as a thin layer of thickness t = 0.5 nm and 
characterized by a dielectric function 1 -i- ia/{cot£o), where a is the surface conductivity of 
graphene, £o is the vacuum permittivity, and CO is the light frequency. Note that the simulation 
results do not rely on the thickness of the film as long as the meshing is sufficiently fine. In this 
paper we consider the case where the photon energy is always less than 2Ef, and under this cir- 
cumstance the intraband contribution to the conductivity dominates. The intraband conductivity 
CTintra of thc graphcnc film can be taken from the standard random-phase approximation |[T9l |20l 



as follows 
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where RbT is the thermal energy, and T is the relaxation time. T can be evaluated by T = 
flEf/eVp, where Vf ^ c/300 is the Fermi velocity and /i = 10,000 cm^V~^s~^ is a typi- 
cally measured, impurity-limited DC mobility [[llO. Then T is approximately 3 x 10~^^ s for 
Ef =0.3 eV. In this work, we consider Fermi energy of 0.3 eV for all calculations, and the 
mode calculations have been performed by using a commercial available finite-element-method 
software package (COMSOL). 

Figure [T] (a) shows the effective refractive indices ^eff of the GPP supported by freestanding 
nanoribbons for two different widths (d = 20, 30 nm). In the frequency range of our interest, 
the nanoribbons with the widths of 20 nm and 30 nm only support the fundamental waveguide 
mode. It has been shown that a wide graphene ribbon may support many modes including edge 
modes and "bulk" modes (TtKTsI. With the decrease of the ribbon width, there exists cut-offs 
for the high-order modes. Eventually, the nanoribbon only has the fundamental mode with even 
symmetry, similar to the dispersion of the surface-plasmon in a metal strip [21]. The effective 
refractive indices shown in Fig. [I](a) of the fundamental GPP of the nanoribbons have large 
values, e.g. 48 for the case of Aq = 15 /im, and d = 30 nm, indicating that the GPP mode 
is extremely confined near the nanoribbon. This is also confirmed by the mode profile shown 
in the inset of Fig.[T](b). The effective refractive index becomes larger when decreasing the 
nanoribbon width from 30 nm to 20 nm. Compared to the nobel metals, graphene generally has 
a relatively large conductivity, thus leading to a lower propagation loss. Propagation lengths L 
as a function of the wavelength are plotted in Fig.[T](b). As an example, when Aq = 10 /im, and 
d = 30 nm, L reaches 0.12Ao, equivalent to 1.2 /im. The propagation loss of the GPP becomes 
weaker when the wavelength decreases, which is mainly attributed to the optical property of 
the graphene material. On the other hand, the propagation length becomes shorter by shrinking 
the nanoribbon width for a specific working wavelength. It can be well understood by the fact 
that when the effective refractive index of the GPP mode becomes larger, there will be a larger 
field overlap with the graphene, thus leading to a higher propagation loss. Similar to surface 
plasmons supported by a metallic waveguide, there is also a competition in the graphene ribbon 
waveguide between the mode confinement and propagation length ||22l [23] |24l |25l . 



3. Bends in graphene nanoribbon waveguides 

We now turn to the question whether waveguide bends (breaking translational invariance) lead 
to reflection or excess absorption? As mentioned above, the graphene nanoribbon waveguide is 
generally lossy. In order to evaluate the performance of the bends, we calculate the transmission 
coefficient and normalize it to that of a straight graphene ribbon waveguide with the same 
propagation length. Propagation properties of the graphene nanoribbon waveguide are analyzed 
by a finite-integral-method software package (CST microwave studio). In all cases, the width 
d of the graphene nanoribbons is chosen to be much smaller than the wavelength, so that only 
the fundamental waveguide mode exists, see the discussion above. 

Figure [2] presents the normalized transmission through a 90^ sharp graphene nanoribbon 
bend, therefore the additional loss associated with reflection and excess scattering loss for the 
sharp graphene waveguide can be obtained. One can see in Fig. [2] that there is no additional 
loss at longer wavelengths. Note that a similar high transmission through sharp bends can also 
be realized in the photonic crystal systems ||26l [27] |28l . For the case of d = 30 nm, the limiting 
wavelength at which the transmission coefficient decreases below 99% is 14 /im. The limiting 
wavelength is blue-shifted to 12 /im when decreasing the nanoribbon width from 30 nm to 
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Fig. 2. Transmission spectra of graphene nanoribbon waveguide bends for d = 20, 30 nm. 
Inset illustrates its top-view corresponding structure and the top-view magnetic amplitude 
at _y = when Xq = 15 jUm and d = 30 nm. 

20 nm. To evaluate the transmittance through the sharp graphene ribbon, we can also employ 
the analogy between the subwavelength waveguiding in photonics and electronics and use a 
common quasistatic approximation ||29j |30l . With this analogy and quasistatic approximation 
(valid as long as d <C A), we can consider the sharp bend equivalent to a junction between 
two transmission lines with the same characteristic impedance. Within the quasistatic approx- 
imation no additional loss emerges for the bend. It is naturally understood that the operating 
wavelength range for having additional bend loss widens with the decrease of d. The inset in 
Fig. m shows the top- view magnetic field amplitude distribution at j = when light passes 
through the sharp comer for the case of d = 30 nm at Aq = 15 /im. To illustrate this behavior 
clearly, we neglect the material loss of the graphene. One can see that light can be perfectly 
turned by 90^ without any noticeable back- scattering or radiation. 
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Fig. 3. Transmission spectra of graphene nanoribbon T- shape splitters (shown in the inset 
with the top view) for d = 20, 30 nm. In this case we choose din = dout = d. 



4. Splitters in graphene nanoribbon waveguides 

Next, we focus on the propagation property of graphene nanoribbon spHtters, which is another 
key element. The calculated normalized transmission coefficients as a function of wavelength 
for graphene T-shape nanoribbon waveguides are plotted in Fig. [3] for the cases of d = 20, 
30 nm. The inset shows a top-view of the structure. Similar to the result for the nanoribbon bend, 
the transmission at the output graphene waveguides becomes larger as the wavelength increases, 
and it saturates at the longer wavelengths, reaching 42%. We find that the limiting wavelengths 
are almost the same as those for the graphene bends shown in Fig.[2l Note that the widths for 
the input and output graphene waveguides have the same values as d. Due to the symmetry of 
the structure, the transmitted power is equally distributed between the two output waveguide 
branches, so that the total transmission is close to 84%. With the quasistatic approximation, 
we can replace the splitter configuration by an equivalent network consisting of a junction 
of three transmission lines with the same characteristic impedance Zq. The combination of 
the two output transmission lines are connected to the input transmission line, thus resulting 
in the equivalent load impedance with a value of = 2Zo. The transmission coefficient at 
each output governed by 1 — [(Z^ — Zo)/(Zl + Zo)]^ is close to 88%. The result obtained by 
the transmission lines has a small deviation when compared with our simulated result. When 
using the transmission line theory, one has to assume that the energy of the guided mode is 
mainly confined in the waveguide and that the distribution of the transverse electromagnetic 
field is nearly uniform within the waveguide. This is only qualitatively correct for the GPP 
supported by the graphene nanoribbon waveguide. However, we can still use the transmission 
line theory to qualitatively interpret the behavior observed in our simulation. In order to improve 
the transmission coefficient of the graphene nanoribbon splitter, we can adjust the characteristic 
impedance of the input waveguide Z/„ by changing the width dj^ of the input waveguide. Figured 
illustrates the calculated reflection coefficient of the graphene T-shaped splitter at Aq = 12 jim 
as a function of djd^^,, where dout, the thickness of the output waveguide branches, is fixed at 
20 nm. An optimal structure with nearly-zero reflection coefficient is observed at djd^^, = 2.25. 
The lower-inset in Fig.|4]shows the top-view magnetic field amplitude distribution at j = when 
light passes through the T-shape configuration for the case of din/dout = 2.25 at Aq = 12 jim. 
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Fig. 4. Reflection coefficient of a graphene nanoribbon T-shape splitter as a function of 
din/dout at A = 12 /im. Inset shows its top-view corresponding structure and illustrates zero- 
reflection when din/dout = 2.25 at A = 12 /im. 



5. Summary and discussion 

We have investigated the performance of bends and spHtters in graphene nanoribbon waveg- 
uides, where the nanoribbon width is much smaller than the working wavelength. We have 
shown that sharp bends and splitters with no additional loss can be designed for longer wave- 
lengths. The effective characteristic impendence model has been used to qualitatively interpret 
the simulation results. As discussed above, the GPP in the graphene nanoribbon has an ex- 
treme confinement, therefore making it difficult to couple light into the nanoribbon by use of 
the butt-coupling technology. Similar to the surface plasmon supported by the noble metals, a 
grating coupler could be used to excite the GPP in the nanoribbons. Recently, the propagating 
and localized graphene-plasmon polaritons have been visualized in real space by scattering- 
type scanning near-field optical microscopy ||3T][32l. Our results further advance the potential 
of such devices by providing building blocks to construct graphene circuits for future ultra- 
compact terahertz devices. 
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